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ABSTRACT 


The purpose of this work Is to highlight the particular aspects of 
miniature heat exchangers performance end to determine en appropriate 
design approach. A thermodynamic analysis Is performed to express the 
generated entropy as a function of material and geometric 
character i st lcs of tho hoat exchangers. This expression Is then used to 
size miniature heat exchangers. 

Key words: Miniature heat exchanger design 5 laminar flow heat 
exchengcr optimization; second law analysis. 


1 . Introduction 


In recent years a large number of applications have been developed for sstall 
superconducting devices being cooled around 4.2 K . Thslr viability will depend on the 
availability of closed cycle, efficient refrigerators. An extended survey CO has 
sho^n that small scalo refr Igerators, for loads around 1 watt, are not cossserc I e H y 
aval I able. 

Different approaches have taen undertaken to build reliable end efficient closod- 
cyclo refr I gsrator 5 . Ih particular, Vul Hosier, G! f fcrd^cf4afcon and Stirling cyclos 
were I nvest Igsted. These cycles are efficient at temperatures greater than 10 K. Tha 
main problems In these prossure-cycl Ing refrigerators for liquid helium temperatures 
are rooted In the poor heat capacity of regenerators materials at low temperatures. 
Preliminary studies C?J snow that for small loads Collins cycle refrigerators are 
advantageous. 

As part of a program to develop an efficient, closed cycle Collins type 
refrigerator for l waft at 4.2 K, we Investigate the performance of miniature helium 
plate heat exchangers (Fig. 1). Due to the rasss flows Involved, these heat exchangers 
operate In laminar f low conditions. 

We will show that the axial conduction In the solid wall has a considerable Impact 
on the performance of low temperature miniature hoat exchangers. 


212 




2. Thermodynamic analysis 


The three major confr ibut Ions to Irreversible it les In heat exchangers are: 

2.1. Entropy generation due to friction, S 

2.2. Entropy generation due to temperature^dl fference 
between the heat exchanger streams, S_ 

2„3. Entropy generation due to heat conduction In the 
solid wall, $ 

c 

Their sum will furnish the generated entropy for the isolated heat exchanger: 


S = S + S T + S . u > 

, . P T c 

2.1. Entropy generation due to friction 

Assuming that helium behaves as a perfect gas, we find that the entropy generation 
due to flow friction Is given by: 

S p = mR[(ip L /p L ) + (AP H /P H >3 , (2) 


where: m = mass flow 

R » gas constant 
p = low pressure 
Pu h *9 h pressure 


Since the high pressure stream friction loss is negligible, S 


P 


is given by: 


S p = "* r V P L } ’ 


(3) 
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where: 


A p. = 0 2 /< «B 2 0 2 )](f/2Re)(U/0) 
f/4Re= frlctlbn factor C3j 
Re = Reynolds number 

p « density 

B,D,L C geometric parameters (Ffg._X). 

Substituting Ap^ Into eq.. (3), we find: 

s p = (m 2 R/p L )(f/2Re)CL/(B 2 0 2 )3 , (4) 

It may be Immediately seen from eq. (4) that entropy generation by friction will 
Increase with Increasing heat exchanger length L and decreasing breadth B and wal l-to- 
wal I d l stance 0, 


2.2. Entropy generation due to finite temperature difference 


The generated entropy due to finite temperature difference Is: 


S T = m 2 c p 2 L(2D/(k f Nu) + f/k^LtTpT^/fTjTylin/fBU) .... 


(5) 


where; c = specific heat at constant pressure of the fluid 
= fluid thermal conductivity 
k r = wall thermal conductivity 
Mu * Mussel t number. 


The derivation of eq. (5) is summarized In the Appendix. 

It may be seen from eq. (5) that the entropy generation will decrease with 
decreasing wall-to-v/al! distance D and wall thickness t and Increasing plate breadth B 
and heat exchanger length L. 


2.3. Entropy generation due to heat conduction In the solid wall 


The temperature gradient In the heat exchanger walls Is linear and the heat 
conduct! vity. k and the crossect ional area A are constant; then the entropy generation 
S associated with the heat transfer by conduction, Q Is: 


s c - WV^iV * (6) 

where: 


’c'.^rVW • (7) 

It may be seer from eqs. (6) and 17) that the entropy generation by axial wall 
conduction will Increase with Increasing crossect lonal area A,*3Bt and decreasing heat 
exchanger length L. * 


2.4. Overall entropy generation 


Tho generated entropy In a laminar counterflow heat exchanger, assuming helium 
behaves as perfect gas end with constant material properties, Is derived front eqs. 
(1), (31,(4), <5) and (6): 
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S => m 3 RCf/(2^Rep L )3CL/(B 2 0 5 ;] + 

m 2c p 2 C<20/k f Nu) ♦ + (8) 

k w*' <T l" T 2 ,2/<T | T 2 ) ^ (Bt/L) ‘ 


It will be useful, for the purposes of simplifying the expression of S end to 
recognize f amlUar dimensionless parameters, to obtain a dimensionless form for eq. 
(8), This Is achieved by dividing eq. (8) by [Ik*/ ^3, that Is constant at pressures 
below 20 atm [43* Equation (8) then becomes: r 


S/C<k f /u>mj = CO' w /k f )it/L)(l.5/Re)tT l -T 2 ) 2 /a i T 2 ))J +■ 

C» (T. -T, ) 2 /(T,T_) ) (Pr 2 Re/Nu) (D/i- >3 + 

1 * * z (9) 

L(k f /k t) )(t/L)(RePr 2 /2)((T ) -T 2 ) 2 /(T,T 2 ))3 + 

C(R/c p )(PrReXf i?/(6p L P L >Xt/D 5 )3 . 

This is the entropy S generated In a counterflow laminar heat exchanger, valid for 
ideal gases In the laminar flow range with Re larger than 100 L3]. 


3. Heat exchanger sizing 


The geometrical parameters chosen aro the wall thickness t, the breadth B, the 
wall-to-wall distance 0 and the heat exchanger length L. Material properties are 

defined by k . 

w 

Results will show that the entropy generation decreases with decreasing t and 

Increasing L and 3, with k 21 k IB) and D « D(L). •“ 

~ a w 

3.1 General methodology for miniature heat exchanger sizing 


Equation (9) shows that tho wall thickness t needs to be small for minimum entropy 
generation. Consequently Its value will be determined by structural requirements. 

The first derivative of S with respect to U Is set to zero: 

w 


« 0.32RePrk f . (10) 

Equation (10) has a form that (s similar to the optimized wall thermal 
conductivity for a concentric tube heat exchanger C53« 

For B considerably larger than D, the Reynolds number may be expressed as: 

Re * 2m/(8u> . (11) 

Substituting eq. U0) Into eq. (9), we find: 

S/Dk f X> >«3 = 1.22Pr((T ) -T 2 ) 2 /(T 1 T 2 ))(t/L) + 

(Pr Z Re/Nu)(T 1 -T 2 ) 2 /(T 1 T 2 K0/L) + (12) 

(PrRe)(R/c p )(fu 2 /(6p L ^))(L/0 3 ) . 
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The minimum of the entropy generation with respect to the wall-to-wall distance 0 
Is then determined slml larly: 


D “ CUR/CpH* » /(2p£ L ))(Nu/Pr))/((T l -T 2 ) 2 /(T 1 T 2 ))3 25 CL3- 5 . (13) 


Equation (13) Is substituted Into eq. (12); S Is then expressed as e function of L 
end Re, where Re Is defined by eq.(ll): 


S/C(k f /> )m] = Cjtt/L) + c 2 Re(ir * 5 , 


(14) 


where: 

c, = 1.22Pr((T 1 -T 2 ) 2 /(T J T 2 )) 

c 2 » 1.05(Pr’ ' 75 /Nu ,75 X(T ] -T 2 ) 2 /(T ) T 2 )) ,75 (R/c p )’ 25 (fu 2 /(p L o L ))* 25 . 

Tho function S Is monotonic decreasing with Increasing L and decreasing Re. The 
remaining parameters B and L may be evaluated by considering space limitation, heat 
leak to tho surrounding and manuf actur I ng costs. 


3.2. Applications for helium heat exchangers 


tfe take Into consideration a laminar counter flow plate host exchanger for a two- 
oxpendGr-saturatod-vapor-pressure cycle refrigerator C6] with a cooling load of 1 watt * 
at 4.2 K. With a reasonable expander efficiency, the mass flow In the balanced heat 
exchanger between 300 K and 60 K Is 0.1 g/s, while tho mass flow In the heat exchanger 
between 30 K and 14 < is 0.035 g/s. 

The generated entropy as a function of the length L, using the Reynolds number as 
a parameter, Is shown In Figures 2 and 3 for the two boat exchangers. 

The heat exchangers breadth B Is computed from eq. (10) for two Reynolds numbers : 
Figure 2: for Re » 100, B_® 0.171 m and for Re « 300, B ® 0.057 m 
. FJgure 3: for Re > 150, B. » 0.294 m and for Re c 300, B e 0.147 m. 

The optimized wall-to-wall distance 0 is computed from oq. (13): 

Figure 2: for G.tm< L < t .Om 0.26nm < D < 0.82mm 

Figure 3: for 0.1m< L < 1.0m 0.06mm < D < 0.19mm . 

Figures 2 and 3 also show the entropy generation duo to axial wall conduction. Its 

contribution to the total entropy generation Is particularly significant for the cold 

miniature heat exchanger and. In general, for low Reynolds numbers. 
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0.0 0.2 0.4 0.6 0.0 1.0 

heat exchanger length L [m] 


Figure 2: Entropy generation rate vs. 

for tfc’o Reynolds numbers 
<T * 300 K, T * 60 K, p. 
Cm « 0 . \ g/s, r * 0.2 rai) L 


length of heat exchanger 
® 2 bar, p u « 20 bar) 

M 



heat exchanger length L [fl] 


Figure 3; Entropy generation rate vs. length of heat exchanger 
parameter: Reynolds number Re 
(T » 30 K, T * 14 K, p =2 bar, p = 20 bar) 

Cm « 0.005 g/S, t « 0.2 fem) H 
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4, Conclusions 


The thermodynamic analysis performed In this work consists ofs 

- expressing the entropy generation contributions due to the axial wall conduction, due 
to the pressure loss and due to the finite temperature difference ss a function of four 
geometrical parameters and one material characteristic. 

- minimizing the generated entropy and thereby reducing the set of Independent 
parameters. 

In miniature heat exchangers, axial conduction In the wall may became an Important 
phenomenon. 

It is shown that the miniature heat exchanger design may not be performed by 
simply scaling down standard heat exchangers because of the short heat exchanger length 
and the laminar flow. 
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Appendix 


The temperature difference AT between the streams compared to the axial 
temperature difference (T.-T^) Is assumed to be small, then the heat transfer rate 
through the heat exchanger surface A Is: 


Q « me (T-T 0 ) » hiTA . (A1 ) 

P * 2 

Assuming tho same convective heat transfer coefficients h for the high and low 
pressure streams, we find: 


h = (NukJ/0 (A2) 

co i 

Thus tho overall heat transfer coefficient h Is: 


t/h = 2/h + t/k . (A3) 

co w 

Considering the fluid heat conductivity k Independent of the pressure, eq. (A2) 
shows that the wall-to-wal I distance 0 Is equal for both channels. 

The finite temperature difference Is derived from eqs. (A2),(A3) and (AD: 


AT = (me /A)<T,-T,,)[<2D/(k # Nu)) + (t/k )3 . (A4) 

p I 2 r w 

Tho entropy generation due to the finite temperature difference AT Is: 

S T = me (ATIT -T 0 )/(T,T-> . (A5) 

1 p l l 12 

By substituting eq. (A4) Into eq. (A5), tho entropy generation becomes: 
S T = m 2 c p 2 C(2D/(k f Nu)> + (t/k^XiTj-T^/lTjiyjD/lBL)] . (A6> 
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Nomenclature 


A area 

B plate width 

Cj constant (n eq. (14) 


c 2 constant In eq. (14) 

0 ^ specific heat of gas 


Re Reynolds number 

S entropy 

t plate thickness 

T temperature 

A difference 

v viscosity 

p density 


D 

f 

h 

h 

c< 

k 

L 

m 

Nu 

P 

Pr 

0 

R 


plate spacing 

friction factor 

heat transfer coefficient 

convective heat transfer coefficient 

thermal conductivity 
heat exchanger length 
mass 

Nussett number 
pressure 
Prandtl number 
heat transfer rate 
gas constant 


Subscr I pts 

1 warm end of exchanger 

2 cold end of exchanger 

c wall conduction 

f fluid 

H high pressure passage 

L low pressure- passage 

P frictional pressure drop 

s cross section 

T gas heat transfer 

w wall 
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